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We calculate the two-body decay rates of "quirkonium" states formed from quirks that acquire 
mass solely through electroweak symmetry breaking. We consider SU(N)i c infracolor with two 
flavors of quirks transforming under the electroweak group (but not QCD) of the Standard Model. 
In one case, the quirks are in a chiral representation of the electroweak group, while in the other 
case, a vector-like representation. The differences in the dominant decay channels between "chiral 
quirkonia" versus "vector-like quirkonia" are striking. Several chiral quirkonia states can decay into 
the unique two-body resonance channels WH, ZH, it, tb/bi, and fH, which never dominate for 
vector-like quirkonia. Additionally, the channels WW, WZ, ZZ, and W"f, are shared among both 
chiral and vector-like quirkonia. Resonances of dileptons or light quarks (dijets) can dominate for 
some vector-like quirkonia states throughout their mass range, while these modes never dominate for 
chiral quirkonia unless the decays into pairs of gauge or Higgs bosons are kinematically forbidden. 



I. INTRODUCTION 

Quirks are fermions transforming under the SM gauge 
group along with a new strongly-coupled "infracolor" 
group SU(N) ic pp. (Earlier ideas were also considered in 
Ref. [2].) The scale of infracolor confinement, Aj c , is as- 
sumed to be much smaller than the masses of all quirks, 
and so the infracolor-strings have an exponentially sup- 
pressed rate to break. Quirks pairs produced in a collider 
remain in a bound state even when produced with large 
kinetic energies. This leads to several interesting collider 
physics and dark matter applications [Tl I5HT^]. (Other 
work on hidden valley models can be found in [T5HT7] .) 
Certain kinds of quirks have already been searched for at 
the Tevatron by the DO collaboration [18] . 

In this paper, we are mainly interested in quirks 
that acquire mass through electroweak symmetry break- 
ing. This is unlike the original model, Ref. PP, where 
quirks acquired "vector-like" masses independently of 
electroweak symmetry breaking. We are motivated in 
part by the discovery that chiral quirks bound in quirky 
baryons can lead to a viable asymmetric dark matter 
candidate [TO]- We do not, however, restrict ourselves 
to the specific theory or detailed parameter choices of 
[TO] . Instead, we consider general SU(N) ic , and calcu- 
late the meson decay rates for both chiral quirks as well 
as vector-like quirks, demonstrating the experimentally 
distinguishable signatures. 

At this point we should emphasize that only some as- 
pects of quirky physics can be calculated (or simulated) 
with standard collider tools. In general, quirks can be 
produced in a standard collider physics process (for us, 
weak production), but then the pt of the quirks must 
be shed before the quirks settle down into a low-angular- 
momentum state. This "spin-down" process is in general 
non-pcrturbative, with the resulting radiation dependent 
on the relative strengths of infracolor and other couplings 
of the quirks. After spin-down and energy loss, the con- 
stituent quirks annihilate, causing quirky mesons to de- 
cay. It is solely this last step that is our interest in this 



paper. 

The annihilation rate of quirky mesons is proportional 
to the lowest non-vanishing radial derivative of the me- 
son wavefunction at zero relative quirk displacement. 
This is entirely analogous to positronium and quarko- 
nium 19J. For an S state, this is ^(O)! 2 , while for a 
P state, |?//(0)| 2 . At high orbital angular momentum 
L, this wavefunction factor is suppressed. Ref. [TJ esti- 
mated the suppression factor in the annihilation proba- 
bility scaling as (f3/L) L+1 / L, where f3 is the quirk rela- 
tive velocity and L > 0. Therefore, instead of annihilat- 
ing immediately, the quirky bound states are expected 
to emit soft radiation to shed their angular momentumj]] 
As the quirky bound state reaches a low angular mo- 
mentum state (L ~ 1), the constituent quirks ultimately 
annihilate; some quirky meson decay rates for certain 
vector- like quirks have been discussed in [JJ [3] IH [H] . 

This paper is organized as follows. We will describe 
our quirk model in Sec. [IT] Next, we present a qualitative 
understanding of the parametric dependencies of the var- 
ious decay channels in Sec. |III[ We present the formal- 
ism to calculate the decay amplitudes in Appendix [A] 
along with the extensive analytical results for all of our 
quirkonia decay rates in Appendices [5] and [C] Much of 
our results for neutral quirkonia can be obtained from 
earlier results on heavy quarkonia [19] . which we have 
compared extensively. Then, we present numerical eval- 
uations of our results, and discuss their implications, in 



Sec. IV Next, we calculate the decay rates for vector- 
like quirkonia in Sec. [V] comparing and contrasting to 
the chiral quirkonia decay results. We summarize and 
provide a clear explanation of which modes have "chiral 
enhancement" in Sec. IVII We conclude with a discus- 
sion, identifying the major signals that distinguish chiral 
quirkonia from vector-like quirkonia in Sec. |VII] 



The radiation may be in the form of soft photons [5] that can be 
detected as rings in the n — 4> plane in colliders. 
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TABLE I: Quirk quantum numbers. 



II. MODEL AND SETUP 

The model we consider is SU(N)i C with two flavors 
in the representations given in Table [I] This is the gen- 
eralization of the model of Ref. [TUj to N infracolors. 
We assume Ai C <C itiq, and neglect the infracolor con- 
finement contribution to the quirky meson masses. The 
Lagrangian that gives mass to the quirks is simply 

C = XuQHu c + X d QHU c . (1) 

Despite the abuse of notation (Q, u c , d c ), we empha- 
size that our quirks are color singlets. After electroweak 
symmetry breaking, the quirks acquire masses M(j,d = 
\u,dv. Writing the electroweak doublet as Q — (u,d), 
we can write the quirks in terms of four-component Dirac 
spinors U, D 

"-{«)"-{*) (2) 

where U, D have electric charge q — ±1/2. The quirky 
mesons formed from these objects include 

(UU), (DD) neutral mesons (3) 
(UD),(DU) q — ±1 charged mesons . (4) 

There are two interesting regions of parameter space 
satisfying the requirement Ai C <C Mu,d- One occurs 
when one quirk mass is much heavier than the other, 
Mjj 3> Md or Mjj ^> Mjj, such that there is one set of 
heavy neutral mesons, one set of intermediate-mass elec- 
trically charged mesons, and one of set of light mesons. 
In this regime, the heavier mesons generically weak de- 
cay to the lightest mesons (microscopically the heavier 
quirks are weak decaying into the lighter quirks) before 
the quirks themselves have time to annihilate. In this 
regime, the relevant annihilation channels consist solely 
of the lightest neutral mesons. 

The second regime, and the main focus of this paper, 
is when Mjj ~ Md. When the two flavors of quirks are 
very nearly degenerate in mass, all of the mesons given in 
Eqs. ([3|,(|4]) are virtually stable against weak decay. All 
of the quirk pairs within the mesons therefore annihi- 
late well before the kinematically-suppressed weak decay 
occurs. This leads to four distinct "towers" of mesons: 
two sets of neutral mesons and two sets of (oppositely) 
charged mesons. 

The neutral mesons (UU) and (DD) can mix with each 
other through infragluon box diagrams that are superfi- 
cially similar to the Vl^-box diagrams within the SM that 



lead to mixing among the neutral mesons of QCD. How- 
ever, unlike QCD, all of the quirks are heavy, while the 
gauge bosons being exchanged in the box diagram are 
massless. This small mixing is an interesting effect for 
further study. Our meson decay rates are invariant un- 
der U f-> D, and we simply compute (QQ) as if it were 
an exact (UU) or (DD) eigenstate. In practice, there 
may be either a small admixture between these states, in 
which case the mixing angle cancels out in our branching 
ratio calculations, or otherwise for maximal mixing, we 
treat (QQ) as the [(UU) + (DD)]/y/2 eigenstate. 

The quirkonium bound state confining potential in the 
Coulombic approximation is |10) 

V(r) = -- t , (5) 

where a contains the relevant couplings for the quirks in 
our model. When infracolor dominates, this is given by 
a ~ a lc = C 2 (N)a lc = (iV 2 - l)/(2N)a lc . The decay 
widths are proportional to the meson wavefunction when 
the two constituent quirks overlap. The wavefunction 
factors that appear in the decay widths, for S and P 
states are 

\Rs(0)\ 2 =4(^a lc M^ (6) 

\R' P (0)\ 2 = ^(ja ic M^ , (7) 

where M is the mass of the meson of the appropriate 
quirkonia state. 

Implicit in evaluating the wavef unctions, Eqs. ([6]),([7j), 
we have assumed the binding energy is dominated by 
the contributions from the infracolor interaction. This is 
not assumed by our analytic results, which are written in 
terms of the radial wavefunction at the origin. Moreover, 
since our numerical results are concerned with ratios of 
decay rates, the dependence on the wavefunction com- 
pletely drops out of the S state quirkonia decay rates, 
and is not particularly sensitive for P states, as we will 
see. 

Transitions between principal quantum numbers can 
occur, just as in bound state problems of QED. The 
transition between the n = 2 P states to the n = 1 
S states is given by the Lyman-alpha electromagnetic 
transition rate. This was estimated for neutral quirkonia 
to be [ID] 

- ^e 2 3 a em £;iJ(0|r|l)| 2 = \(J^j a em af c M, 

(8) 

where M is roughly the meson mass. Charged quirko- 
nia have the same rate, so long as infracolor dominates 
the potential, Eq. As we will see, the Lyman-alpha 
transition is typically faster than the annihilation rates 
of P state quirkonia. There are exceptions, however, for 
chiral quirkonia, which we calculate below, and discuss 
the resulting final state signatures. 



2 



III. QUALITATIVE DESCRIPTION OF 
RESULTS 

The formalism we used for our quirkonia decay calcu- 
lations is given in Appendix[XJ and the complete analytic 
results are given in Appendices [E] and [C] Here, we eval- 
uate the parametric scaling of the various transition and 
decay rates. There are five qualitatively distinct rates 
involved in quirkonia decay: 



r 








(9) 




r(S -+ g'g') r 






(10) 


T(S- 


SM + SM) r 


- Osm«L m ( 


y TO 2 / 


(11) 




r(P g'g') - 


^aJ c M 




(12) 


r(p - 


-> SM + SM) - 


- aluatM ( 


'Af 2 ^ 

> TO 2 / 


(13) 



The first rate corresponds to the Lyman-alpha electro- 
magnetic transition given by Eq. ^ above. The second 
rate, S — > g'g' refers to specifically the 1 5'o state decaying 
into a pair of infragluons. Note that the 3 Si quirkonia 
state does not decay into a pair of massless gauge bosons, 
due to angular momentum conservation, just like the 3 Si 
quarkonia state [19]. The third rate, S ->• SM + SM 
refers to 1 Sq or 3 Si state decaying into a pair of SM 
particles with SM coupling agM- The fourth and fifth 
rates, P g'g' and P — > SM + SM, refer to any of 
the P states decaying into the above modes. Generally, 
if the Lyman-alpha transition is possible (i.e., for any 
of the P states) , it dominates over the quirkonia decay 
modes. The exception to this is if there is longitudinal 
or Yukawa "chiral enhancement" , which can occur either 
singly {j3 = 1) or doubly (f3 — 2) depending on the final 
state. Double-longitudinal enhancement tends to over- 
come Lyman-alpha emission when M is somewhat larger 
than one of the SM bosons, m = Mw, Mz, compensat- 
ing for the suppression by the larger number of couplings. 
Calculating exactly which modes are enhanced, and why, 
is the main thrust of the paper. 

The residual dependence of the decay branching ra- 
tios on the parameters of the theory, a^ Ai Cl and the 
infracolor group SU(N)i C , arise from: a) neutral quirko- 
nium decay into infraglueballs, b) possible Lyman-alpha 
infraglueball transition, and c) Lyman-alpha electromag- 
netic transition. 

The first issue, neutral quirkonium decay into in- 
fraglueballs, is handled by choosing to evaluate "width 
ratios" into visible SM particles, rather than the stan- 
dard branching ratios of neutral quirkonium decays. This 
is because the infraglueballs are expected to be very long- 
lived and escape the detector [I], and thus annihilation 
to g'g' is expected to yield no hard SM resonance signal. 
The possibility of 3-body decays, going into a pair of in- 
fraglueballs as well as a SM particle, is beyond the scope 
of this paper. 

The second issue, Lyman-alpha infraglueball transi- 



tion, can occur if the kinematics of the transition per- 
mit it. The energy difference between the n = 2 P state 
and the n = 1 S state is (3/32)a 2 c M. The infraglue- 
ball mass is of order, but somewhat larger than A^ c . For 
a ic ~ O(O.l) and A ic ~ 0(1 GeV), the infraglueball 
mass is already close to this energy splitting, and so kine- 
matic suppression is generic for somewhat larger values 
of ai c . For smaller values of Ki C , presumably accompa- 
nied by smaller values of a.i C , the infraglueball emission 
is less kinematically suppressed but has an overall tran- 
sition rate that is smaller. In this region of small ai c , the 
P states are more likely to transition to S states before 
annihilating. A different choice of the infracolor group 
introduces an order one change to <5i C , as well as the run- 
ning of 6ti c in the Coulombic potential |10j . Since these 
changes appear as shifts in the coupling of the Coulom- 
bic potential, they can be at least partially absorbed by 
a redefinition of cn c . 

Finally, there is the different dependence on 6n c (and 
c*Sm) between the Lyman-alpha transition rate versus 
the quirkonia decay rates. The Lyman-alpha transi- 
tion scales as asM&ic whereas the decay rates of the 
P states scale as a| M af c through \R' P (0)\ 2 . Hence, there 
is a relative suppression of branching and width ratios of 
roughly asM&ic relative to the Lyman-alpha electromag- 
netic transition rate. 

Consider two regimes, one in which the Lyman-alpha 
transition rate is dominant and the other where it is sub- 
dominant to decay rates that are doubly-enhanced. In 
the latter regime, the branching and width ratios are 
rather insensitive to oti C , since the only dependence on 
(Xi c enters from the wavefunction. This dependence drops 
out, analogous to the S states. In the first regime, how- 
ever, the total width is dominated by the Lyman-alpha 
transition rate. The width and branching ratios of de- 
cay processes scale linearly with a.i C - But, the regime 
where the potential is dominated by infracolor while still 
allowing a Coulombic approximation to the bound state 
potential only allows for about one order of magnitude 
change in oti c . 



IV. BRANCHING RATIOS AND WIDTH 
RATIOS 

We now present our results for charged and neutral 
quirkonia decay rates. In what follows, we carry out 
several numerical computations of branching ratios and 
width ratios, using the results from Appendices [B] and [C| 
to demonstrate the dominant SM decay channels for the 
various bound states of quirkonia. 

In the following, we have chosen a specific infracolor 
group, N = 2, and infracolor coupling, a ic = 0.2. As 
we discussed above, the additional relative suppression 
factor of the P state decay rates is only one power of cxi C - 
The value a, c = 0.2 tends to maximize the possibility 
that P states can, at larger quirkonium masses, annihi- 
late into SM modes before the Lyman-alpha transition 
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A. Charged quirkonia 

We now discuss the branching ratios of charged quirko- 
nia. The analytic results for the decay rates can be found 
in Appendix [B] Given a final state /, the branching ratio 
for /, BR(QQ -> /), is 



BR(QQ -> /) = 



r(QQ -> /) 
£,r(QQ->/) : 



(14) 



where the sum is over all final states. 

The charged quirkonium case is particularly simple. As 
the system is electrically charged, it cannot decay into 
g'g'. Fig. [T] shows the decay branching ratios of charged 
quirkonium states. For all states, only the WH partial 
width is sensitive to different values of the Higgs mass. 
The plots shown for charged quirkonia here are also ap- 
plicable to any massive bound states that only decay via 
the electroweak SU (2) x U(l) group, with electric charges 
Qu = —Qd = 1/2. Note that we only show the summed 
width over the massless fermions (2 quark pairs, 3 lepton- 
neutrino pairs) , as t he w idths of all massless fermions are 
the same (see App. B4). Also, we only show the decays 
for the UD meson. We have checked that the widths for 
UD decay are the same. 

The branching ratios of different bound states are plot- 
ted in Fig. [T] For the S states, the W Z partial width 
dominates. For P states, radiative transition usually 
dominates. For 3 Pi, the WH width becomes larger than 
the Lyman-alpha transition width when the meson mass 
is larger than > 600 GeV, provided that the meson is 
heavier than the threshold. 



B. Neutral Quirkonia 



where eQ is the quirk electric charge. Setting A = 3 
reproduces the results in [19] . The 3 Si state cannot decay 
into g'g'. Its decay into g'g'g' is given by [3], 

TY 3 <? ^ n'n'n'\- iN2 - l){N2 - A) ^ - ^ I ff <fi\\* 

(16) 

This vanishes for the special case SU(2)i c , since three 
gluons cannot form a infracolor singlet. 

Instead we present our results in terms of a "width 
ratio" 
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(17) 



Also, for reasons of clarity, we do not present the plots for 
the branching ratios when the Higgs mass deviates from 
125 GeV. Unless the final states involve Higgs bosons, a 
larger Higgs boson mass would only push the correspond- 
ing thresholds towards higher meson masses, leaving the 
other width ratios mostly unchanged as in the case of 
charged quirkonia. However, there is a qualitative change 
in the width ratios for the 3 Pq state when the Higgs mass 
is sufficiently large, which will be discussed below. 



2. 1 So and 3 Si 



21) 



The width ratios for the S states are shown in Figs. 2a 
The dominant decay channels are ZH and ti. The 



ZH channel receives double enhancement from the lon- 
gitudinal Z and the Yukawa coupling between the Higgs 
boson and the quirk. Also, the ti channel is enhanced by 
the top mass. 

The results for the 3 Si state can be discussed more pre- 
cisely because the infraglue channel is absent. Even for 
values of M not far away from 2mw, the double longitu- 
dinal WW mode dominates. Because of Bose symmetry, 
the two Z , s cannot be longitudinal simultaneously and 
the ZZ mode is suppressed. 



The results for neutral quirkonia are more complicated 
than their charged counterparts. Not only that there are 
more decay channels, but in some cases the mesons can 
decay into infraglue pairs, g'g' , that hadronize into an 
infraglueball pair <fJ<j>'. 
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'Pi 



The width ratios are shown in Fig. 2d is dominated 
by the Lyman-alpha transition throughout the sub-TeV 
range. All other widths contain a single enhancement 
factor, from either the longitudinal mode or the quirky 
Yukawa. 



As the mesons are color singlets, only t- and w-channels 
contribute to the decay amplitudes of B — > g'g' . Then 
the decay rate should be proportional to that of B — > 77. 
A simple calculation shows that, for an SU(N) ic color 
gauge group, 
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(15) 



4- 3 Po 

The 3 Po width ratios exhibit an interesting feature 
when the Higgs mass is larger than 2mw, 2m z and 2M t , 
where M t is the top mass. The decay channels WW, 
ZZ, and ti involves an s-channel Higgs boson exchange. 
When the meson mass is near the Higgs mass M ~ Mh, 
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FIG. 1: Decay branching ratios of charged chiral quirkonia in different J PC states. Solid lines are with Higgs mass 

M H = 125 GeV, dashed lines with M H = 250 GeV. 



the widths are enhanced by the s-channel Higgs reso- 
nance. This can be seen in Fig. [3] There, the WW and 
ZZ widths have a resonance at M = Mh = 250 GeV 
when the s-channel Higgs boson is on-shell. The it width 
does not exhibit this behavior because at 250 GeV, the 
decay into two top quarks from a single Higgs boson is 
kinematically forbidden. 



-Pi 



The branching ratios for the 3 Pi state are shown in 
Fig. [2e] The ZH channel are doubly enhanced and is 
dominant for M > 700 GeV. 
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FIG. 2: Decay width ratios of neutral chiral quirkonia in different J PC states. Solid lines correspond to a Higgs 
mass Mfj = 125 GeV, while dashed lines correspond to Mh = 250 GeV. In many instances, there is no difference 
between the width ratios for different Higgs masses, and thus the solid lines overlap the invisible dashed lines. For 
figure (c), we have presented the choices Mjj = 125 GeV. We illustrate the difference in decay width ratios changing 

to M H = 250 GeV in Fig. |j 
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COMPARISON TO VECTOR-LIKE 
QUIRKONIA 



Annihilation rates for the case of vector-like quirks in 
The channels WW, ZZ and HH are doubly enhanced certain other representations has been calculated in [4. 
and will take over the radiative transition at high meson ,13.. There is not a general rule that relates the decay 
mass (> 1 TeV). rates of vector-like quirks to chiral quirks. But in certain 
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FIG. 3: 3 P with M H = 250 GeV. The resonance 
structure is due to the s-channel Higgs boson. 



circumstances one can be obtained from the other, and 
vice versa. In this section, we will discuss differences and 
similarities of vector-like and chiral quirk decay rates, 
and give examples in cases where the decay rates are 
related. 

We wish to compare our results for chiral quirks to a 
related theory with vector-like quirks. The vector-like 
theory we consider consists of the doublet Q given before 
in Table |TJ but now we replace 



u c (N, 1,-1/2) 
d c (N,l,+l/2) 



Q'(N,2,0) 



(18) 



Yukawa couplings, Eq. ([I]), are not present, and instead 
we write the vector-like mass MQQ' where M = Mjj = 
Md- There are several differences that lead to qualita- 
tively different decay widths. 

First, the coupling of electroweak gauge bosons to left- 
and right-handed quirks are the same - the quirk- W/Z 
coupling is a purely vector interaction, and processes 
that proceed through the axial vector coupling in the 
chiral case are absent for vector-like quirks. As an ex- 
ample, consider the decay rate r( 3 Si — > //), for neu- 
tral and charged quirkonia. The only diagrams are the 
s-channel j/Z or W . In the neutral case, the only differ- 
ence that separates vector-like and chiral is the different 
axial-vector and vector coupling of the Z. Therefore, 
the expressions for vector-like [4j and chiral [19] are the 
same. For the charged case, the axial-vector and vector 
couplings are not explicitly written in [3] , but their rate is 
4 times larger than the chiral case in [T§] . This is because 
the s-channel W couples to both left and right handed 
quirks in the vector-like case, whereas in the chiral case 
they only couple to left handed quirks. Therefore, the 
decay rate into a fermion-antifermion pair for a charged 
Si is four times larger than its chiral counterpart. 

Second, the quirks do not couple to the Higgs and 
the corresponding Goldstone bosons (through the lon- 
gitudinally polarized electroweak gauge bosons) . Virtual 
Goldstone bosons can only appear in the s-channel, and 
since the Goldstone bosons are pseudoscalars, they only 
contribute to the 1 S$ decay rates. Vector-like quirks, 



by contrast, do not have couplings to the Higgs or the 
Goldstone bosons. In addition, Goldstone bosons can 
appear in the final state (appearing as longitudinally po- 
larized electroweak gauge bosons). This leads to quali- 
tatively different decay rates into gauge bosons for all of 
the bound states. 

For completeness, we present the width and branching 
ratios of vector-like quirkonia in Fig.[4]for neutral quirko- 
nia and Fig. [5] for charged. There are striking differences 
between the chiral and vector-like cases. The most promi- 
nent feature in the vector-like case is that all decay widths 
have the same asymptotic behavior at large quirkonium 
mass - there are no longitudinal enhancements of W/Z 
anywhere. This is expected, as the longitudinal W/Z 
asymptotes to the respective Goldstone bosons, which do 
not couple to the vector-like quirks in u- and i-channcl 
quirk-exchange diagrams. Also, the trilinear gauge boson 
coupling appearing in s-channel gauge boson exchange 
arises from the electroweak gauge structure of SU(2)l 
and has no relation to the electroweak breaking mech- 
anism. Therefore, one would not expect any enhance- 
ments in the decay widths of vector-like quirkonia. With- 
out longitudinal enhancements, the Lyman-alpha tran- 
sition dominates over all P-state decays for all quirko- 
nium masses. Whereas in the chiral case, decay channels 
that receives longitudinal enhancements can dominate 
the Lyman-alpha transition at large quirkonium masses. 

In the low quirkonium mass regime, the overall be- 
havior of both vector-like and chiral quirkonia are simi- 
lar: P-states predominantly decay via the Lyman-alpha 
transition and 3 Si into qq. It is interesting to note that 
for So, jZ dominates the vector-like quirkonium decay, 
whereas 77 is dominant for chiral quirkonia. This is be- 
cause the primordial electroweak gauge boson W£ cou- 
ples not just to the left-handed vector-like quirk, but to 
the right handed one also! A rough estimate indicates 
that this gives a factor of four increase in the ^yZ rate for 
the vector-like case. Indeed, the isospin contribution to 
the vector coupling of the Z to the quirks for the vector- 
like case is twice as much as that for chiral quirks. 



VI. UNDERSTANDING CHIRAL 
ENHANCEMENTS 

Figures [T] and [2] show that decay processes can be singly 
or doubly enhanced by either the Yukawa coupling or the 
longitudinal modes of gauge bosons at high quirkonium 
masses. A summary of enhancements received in the de- 
cay channels can be found in Tables [H] and |III| In this 
limit, the Goldstone equivalence theorem applies and the 
enhancements in various decay channels can be seen by 
the matching of the J, P, and C numbers between the 
decay quirkonium state and the final state consisting of 
Goldstone/Higgs bosons and/or transverse gauge bosons. 
First, we determine the J PC numbers to the final state 
particles; h and ++ to the Goldstone boson and Higgs 
boson, respectively. For photons and transverse W and 
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FIG. 4: Decay width ratios of neutral vector-like quirkonia in different J PC states. Solid lines correspond to a Higgs 
mass Mfj = 125 GeV, while dashed lines correspond to My = 250 GeV. In many instances, there is no difference 
between the width ratios for different Higgs masses, and thus the solid lines overlap the invisible dashed lines. 



Z, we assign P and C numbers according to the P and C 
of the bilinears ipj^tj} and ^pj^^ip. The photon has J PC 
= 1 , and the W/Z has 1 for the vector coupling 
and 1 ++ for the axial vector coupling (that is, we absorb 
the violation of C and P from the axial vector coupling 
into the J PC of the W and Z). For the quirkonia, we 
have CT+ for 1 S , \— for 3 5i, l+~ for 1 P 1 , 0++ for 3 F , 
1++ for 3 P 1 and finally, 2++ for 3 P 2 . 

Next, we write down all the available J PC with differ- 



ent orbital angular momentum L between the two final 
state particles, and match with the J PC of the quirky 
meson to determine which meson decay channels are en- 
hanced. As an example, consider QQ — > Z7, where 
there are only the t- and u-channel diagrams. Single 
enhancement from the longitudinal Z is present in some 
meson states. In the limit where the quirkonium mass 
M 3> Mz, the longitudinal Z is equivalent to the corre- 
sponding Goldstone boson <jp . The J PC of the </>°7 sys- 
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FIG. 5: Same as Fig. [4j but for charged vector-like quirkonia. Only the two J states shown here have nontrivial 
branching ratios. The state 3 Si cannot decay through two-body decays. The 3 -Po,i,3 states can only decay 

radiatively into S states. 
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TABLE II: Decay channels that receive enhancements 
for charged quirkonia. 





Singly enhanced 


Doubly enhanced 


'S 




ZH 


3 Si 


Z-y,ZZ, ZH, jH 


WW 


'Pi 


Zj, WW, ZZ, ZH, -yH 




3 Po 




WW, ZZ, HH 


3 Pi 
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Zj, ZH 


WW, ZZ, HH 



TABLE III: Decay channels that receive enhancements 
for neutral quirkonia. 



tern is determined by combining that of the Goldstone 
boson (0 h ) and that of the photon 1 , which gives 
P . Keeping the total angular momentum J < 2, we 
can add orbital angular momentum L into the system, 
forming {0 -- , l -- , 2~ } for L = 1, and {l + -,2+~} for 
L = 2. This is captured in Table IV One sees that, only 
the states l" 1 and 1 matches with the existing quirky 



meson sta tes 1 P\ and 3 Si, respectively. One can see from 
Appendix C2 that only the 1 P\ and 3 Si states are en- 
hanced. One can also see that both of the decay proceed 
via the axial vector coupling of the Z. This is because 
the 4>° ff coupling is a pseudoscalar coupling ipj 5 ^, and 
this must correspond to the axial vector coupling of the 
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TABLE IV: J PC of the 4>°~/ system. 
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TABLE V: J PC of the Z^J system with vector coupling 
of the Z. 



Z, V7 P 7 5 '0- 

We will go further and illustrate that all other meson 
states decay to Zj via the vector coupling and receive 
no enhancements with the same procedure. Consider the 
final state with a transverse Z and 7, both with J PC 
I for the vector coupling of the Z. The J PC of the fi- 
nal states are collected in Table [V] where one can see that 
the 1 S'o, 3 Po 3 Pi, and the 3 P^ states are not enhanced. 
The J PC of the final states with an axial vector coupling 
can be obtained by flipping the C and P numbers ev- 
erywhere in Table |V| Then one sees that there are also 
contributions with no longitudinal Z enhancements via 
the axial vector coupling to the decay widths of 3 Si and 
1 Pi. Indeed, there are terms in the expressions for the 
corresponding decay widths that are not enhanced. 

The procedure above can explain a large number of 
enhancements for different quirkonium decay processes. 
However, there arc instances where the procedure pre- 
dicts leading enhancements in some processes when there 
should not have been any. The fictitious leading enhance- 
ments predicted by this procedures are: 1 Sq — > W^<p + , 



X -> 0+7, x Pi -> + P, 3 P -> W+H and ^So, 
3 P } -)■ Z T <£+. 

The 1 S'o state are not singly enhanced in two-body 
systems consisting of one transverse gauge boson. This 
can be seen by tracking the spin along the quirky 
fermion line. First, the 1 S'o projector is proportional to 
Ssis 2 ^iSo 2 " 81 ^ 2 ' wnere Cig* 2 is the Clebsh-Gordon co- 
efficients corresponding to the 1 Sq state, which vanishes 
for si = s%. In the systems being considered, there is only 
one vertex that flips the spin - the fcrmion-fcrmion-gauge 
boson vertex. The rest are either scalar (Higgs boson) or 
pseudoscalar (Goldstone boson) couplings that do not in- 
volve a spin flip. Also, one can see that the Dirac spinor 
v has an opposite spin compared to the Dirac spinor u 
in the relation v s = —2sj 5 u- s . Therefore, the state 
cannot decay into any states consisting of only one trans- 
verse gauge boson (with the other outgoing particle being 
the Higgs or Goldstone boson) and the 1 Sq decay width 
cannot be singly enhanced. 



VII. DISCUSSION 

We have calculated the chiral quirkonium decay rates 
for quirks that acquire mass through interactions with 
the Higgs boson. While in this paper we have not stud- 
ied the production rates of quirks at the LHC, this is 
straightforward using standard collider physics tools, and 
has been done recently in the literature. For vector-like 
quirks with clcctrowcak quantum numbers, the produc- 
tion cross section was calculated in Fig. 6 of Ref. [15] . 
There it was shown that pair production of uncolored 
quirks produced through j/Z exchange has a cross sec- 
tion ranging from roughly ~ 10 5 to ~ 10 fb at y/s = 7 
TeV LHC for the quirkonium mass range of 100 to 
1000 GeV. We expect that the chiral quirk production 
cross sections are very similar in size. Given the spectac- 
ular signals that result from quirkonia decay, the discov- 
ery of quirkonium resonances can occur quickly. This is 
likely to occur well before the detailed properties (spin, 
P, C, etc.) of the resonances can be determined. 

The particular decay channels not only can tell us 
about the constituent quirks' quantum numbers, but per- 
haps even more interestingly, how the quirks acquire 
mass. Quirkonia with chiral quirks have longitudinal 
and Yukawa enhancements that are absent or highly sup- 
pressed in quirkonia with vector-like quirks. In this pa- 
per we have demonstrated the striking differences be- 
tween the dominant chiral quirkonium decay channels 
as compared with vector-like quirkonium decay channels. 
This should enable the LHC to easily distinguish whether 
quirks are chiral or vector-like from the observation and 
branching ratios of the dominant decay channels. 

For electrically-charged chiral quirkonia composed of 
quirks with the quantum numbers given in this paper, 
the state 1 5'o decays into WH or tb/tb overwhelmingly 



for quirkonium masses larger than about 250 GeV. Con- 
trast this with vector-like quirkonia, where the 1 Sq with 
the quantum numbers given earlier in the paper, does 
not even have two-body decays. Chiral quirkonia in the 
3 Si state have WZ is the dominant decay channel. For 
vector- like quirkonia, the //, summed over all flavors of 
SM fermions, dominates. We also demonstrated that the 
Lyman-alpha transition is dominant in all of the charged 
quirkonia P-states, except for 3 Pi for chiral quirkonia 
when the "doubly-enhanced" WH decay becomes signif- 
icant for quirkonium masses > 600 GeV. The W*y de- 
cay deserves more discussion. It was shown in [5J that 
the channel is dominant when their squirk and anti- 
squirk pair has low relative velocity, in other words, an 
5-state. For our case, the partial width vanishes due 
to our choice of quantum numbers: Qu = —Qd — 1/2, 
where we found the Wj pa rtial width is proportional to 
2 , see Appendix 



ID 



Bl 



(Qu 

For electrically-neutral chiral quirkonia, the dominant 
decay channels of 1 5'o are ZH and tt. Again, contrast this 
with vector-like quirkonia where WW or ZZ dominates. 
For 3 Si, the WW channel dominates for chiral quirko- 
nia, versus // for vector-like quirkonia. For all P-states, 
Lyman-alpha emission dominates for vector-like quirko- 
nia, whereas for chiral quirkonia there are several decay 
channels that can become significant when the quirko- 
nium mass is large. In particular, HH can dominate for 
3 P , ZH for 3 Pl, and WW for 3 P 2 . 

Indeed, perhaps one of the most interesting decay 
channels that we found is the 3 Po decay into two Higgs 
bosons, which becomes the dominant decay channel for 
quirkonium masses > 500 GeV. This could be striking 
signal at LHC, given that the di-Higgs system would re- 
construct to an invariant mass peak of the 3 Po state. 

Finally, it is tempting to consider applications of our 
results to various existing hints at colliders. For instance, 
the prominence of the electrically-charged chiral quirko- 
nium decay channel, 1 S$ —> WH is suggestive: chiral 
quirkonia with mass M ~ 300 GeV with some minor 
modification of the H decay into jets could easily lead 
to the CDF excess in the Wjj signal [HJ [55] . The cross 
section can be easily adjusted to match the excess, sim- 
ply by enlarging the number of infracolors or flavors of 
quirks. Given the incredible performance of the LHC 
over the past several months, we leave this pursuit to 
future work. 
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Appendix A: Matrix elements of bound state decays 



This section reviews the procedures to evaluate the decay amplitudes of different angular momentum bound states 
following the method in |20j . We work in the non-relativistic limit, where the relative momentum of the constituents, 
|q| -C M, where M is the mass of the meson. We also ignore the contribution to the meson mass from the binding 
potential, i.e., we take M = 2toq, with rriQ the mass of the individual quirks. 

Calculations of the matrix element involving an incoming bound state and an outgoing free state, (X\iT\B), are 
needed to evaluate different bound state decay rates. This is most conveniently done by writing the bound state as a 
superposition of free fermion states with spins (s±, S2) and momenta (pi,P2): 



\B) = \ 2s+ %) = J2< lmss ^)\ lmss 

MS Z 

T fd 3 q 



M 7(2tt) 3 



S1P1S2P2}, 



(Al) 



where tp is the Schrodinger wavefunction of the bound state. In its rest frame, p\ = Q/2+q, and P2 = Q/2 — q, where Q 
is the 4-momentum of the meson, and q is the relative 4-momcntum of quirks. Then, the quantity <X\iT\sip\S2P2>= 
iv S2 {p2)^u Sl {pi) is the usual fermion-antifermion annihilation matrix element into the outgoing state /. Expanding 
the above to the lowest non- vanishing order in q, we found the following decay amplitudes for S and P states, 



M 1 So) = xl^^Rs^Tr^i-Q + M)], (A2) 



A{3Sl) = \lY^M Rs{0)Tr[ ^^ + M)] ' (A3) 



■U'/V) = -i\I^R , P (0)Tr[±e^ 7 H-<8 + M)+^ftjt 7 % (A4) 



A ( 3 Po) = hi ^7-MO) Tr[\jt a ( %| - la ) + M) 3^], (A5) 



4ttM ^ w l 2 V M2 
A( 3 Px) = i^^R' P (0)Tr[2^h 5 - -±- e ^ s Q p ^ alp e s (-($ + M)], and (A6) 

A( 3 P2) = -i\[^R'p{Q)Tr[\jt a e aP lp{-® + M)}, (A7) 

where R(0) and i?'(0) are the meson radial wavefunction and its derivative at the origin, respectively. 

The mass of each meson is distinct, and in principle the M in each of the above expressions should be replaced with 
the mass for that particular bound state. Since we assume cii C is perturbative, the differences between bound state 
energies is parametrically small, of order oq c M. In most instances, one can use our expressions below, substituting 
the proper quirkonia mass for M, and computing the rates. For the branching ratio plots we present below, however, 
this difference is small. 

The quantity = /dq a is the derivative of the matrix element with respect to the relative momentum q. 
The meson polarizations in the rest frame, e M for spin-1, and e MI/ for spin-2, are chosen to be 



e 



M _ 



e£ = (0,0,0,1), 

e "f = £<1M,LS,|2J Z ), (A8) 

M,S Z 

with J z G {—2,-1,0,1,2}. The factor of yN arises from normalizing the meson as an infracolor singlet, exactly 
analogous to what is done with the QCD color factor for quarkonia [T5] . 

The P and C parities of the above angular momentum states are manifest in each of the decay amplitudes above. 
For example, with = (M, 0, 0, 0) in the rest frame, the P and C parities of the bilinear constructed from the 
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projector appearing in A( 1 S , ), ip^ 5 (—(£) + M)-f°tp, are — and +, respectively. Thus J PC = CT + for 1 Sp, as expe cted . 
One can check that the other amplitudes give the expected J PC using the same procedure. From Eqs. (A2) to (A7|, 
we rederived all of the two-body decay rates listed in [19 |^] 



Appendix B: Decay Rates of Charged Quirkonia 

Apart from a color factor of N, decay rates of neutral quirkonia that do not involve any gluons are the same as 
listed in [19]. The decay rates of charged quirkonia will be discussed in this section. 

Charged quirkonia are expected to have larger partial widths than their neutral counterparts. This is because 
charged particles do not have a well-defined charge conjugate parity, hence loosening the constraints posed by CP 
conservation. Here, we list the partial widths of charged quirky mesons with positive unit electric charge, i.e. Qu — 
Qd = 1, where Q is the electric charge of either the up-type or down- type quirks. The mass ratio squared Ri and the 
relative velocity Pij appearing in the formulas below are defined as 



= ^Jl + (R l -R J ) 2 + 2(R i + R j ), 
respectively, and cw = cos 9w is the cosine of the Weinberg angle. 



(Bl) 
(B2) 



The charged quirkonium decay widths into W + j are qualitatively different to the widths of neutral quirkonia into 
Z*y. There are two reasons for the differences; the decay into W + j can go through an s-channel with a W exchange. 
The corresponding diagram is absent for Z7; the photon does not couple to a electrically neutral Z. Another reason is 
that the photon couples to quirks of different electric charges in the t- and w-channel diagrams, due to the emission of 
a charged W. It is illuminating to write the chiral projection operators as Pl,r = (v w Taw7 5 )/2, with aw — vw = 1, 
so that the vector and axial-vector contributions from the W current are manifest. The partial widths into W + ^/ are 

Naa w v^ 



TCS+ -> W + j) 



T( 3 S+ 

rfp+ 
r( 3 p+ 
r( 3 p+ 



t( 3 p+ 



4 }[2 (Qu + Qd) 0- ~ R w )\Rs(0)\ , 
Naaw f w { Qu + Q D ) 2 {1 R W )\R S (0)\ 2 , 



12m 2 w 
Naaw 



[a 2 w (Qu + Qd) 2 (1 - R 2 W ) + v w (Qu ~ Qd) 2 Rw(1 ~ R W )]\R P (0)\ 2 , 



Naotwi)- — Rw) 

Naaw 
2Mh 



a w (Qu - QdY + v w (Qu + QdY I 1 



it' 



a 2 w (l - R w ) + 4*%r(Qv + Qd) 2 R 2 w 



1 — Rw 
R' P (0)\ 2 , 



\Rp(0)\ 2 



Naa w (l — Rw) 
10M 2 m^ 



a 2 w (Qu-QD) 2 (3 + 4R w ) 



1 + Rw 
1 — Rw 

4v w (Qu + Q D ) 2 R W {Q + 3R W 



R 



(l-Rw) 2 



Interestingly, all but one term in T( 3 P 1 — > W + "f) are proportional to either the hypercharge Y = (Qu 
the isospin T$u — (Qu — Qd)/2 of the quirks. 



(B3) 
(B4) 
(B5) 

(B6) 
(B7) 

|P^(0)| 2 . 

(B8) 
Qd)/2 or 



2. W + H 

In the limit degenerate quirk masses, their coupling constants to the Higgs boson are the same. As a consequence, 
the decay matrix elements has the same form as that for the decay into ZH , and can be obtained by the replacements 



2 We found a relative sign difference between the two terms in 
the amplitude A( 3 Pi). We attribute this to our definition of 
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gz — > g/V% and vector and axial vector couplings by 1/2; a, v —> 1/2. This gives a conversion factor of 1/(2a/2) 
converting the ZP matrix elements to VFP: 



(B9) 



Therefore, the partial widths into W + H, have exactly the same form as those for ZH, aside from a factor of 1/8. 
The analysis for ZH in [TH] applies to W + H as well. The partial widths are 



r( 1 5 + -> 
r( 3 s+ -> 



r( 1 p+ vr 1 
r( 3 p+ -> W" 1 
r( 3 p+ -> W" 1 



32M 2 



i? 2 



|i?5(0)| 2 , 



(BIO) 



Na 2 w p WH M 2 {8R W [(1 - Rw) 2 + R H {1 - 3R W )] 



m w V 



{l-R w ) 2 {l-R H -R w f 



384 

[P H (1 - 3i?w) — 2i2ff(l - i?i^(2 + + (1 - Piv)(l - Rw — Pwh) 



(1 - R W ) 2 {1 - Rh - Ri 



\Rs(0)\ 2 , 



Na 2 w (3 



WH 



H) AM 2 m 2 w {\ - R H - R w ) 2 
B) = 0, 



H) = 



Na 2 w p WH /2[1 -R H + Rw?[l + Rw(2 - R H + Rw)] 2 



8m^ 



+Rw 



(1-R H -R W ) 2 (1-R W ) 2 

1 2n 



Pwh ~ 4(1 — Rh — Rw) 



r( 3 p+ -> vk+pt) = 



4P W 

1 -P w + (1-R H - Rw) 2 
3Na 2 w P WH 



|P' P (0)| 2 , 



40M 2 m 2 / (l-P ff -P w ,) 4 



IP'p 



(Bll) 
(B12) 
(B13) 

(B14) 
(B15) 



WZ 



Notice that double longitudinal modes are allowed from the decay of a charged quirkonium in the state. This 
is impossible for the neutral quirkonium case, where it decays into ZZ or WW. To see this, the 1 5o state has 
J PC = h , but at zero angular momentum, the double longitudinal state has J PC = ++ . The decay into double 
longitudinal modes for neutral quirkonia in Sq is forbidden by CP conservation. For charged states, the charge 
parity is irrelevant, and the decay into double longitudinal mode is allowed by CP conservation. Naively, one would 
expect the 1 So decay rate is longitudinal from appearance of the \ j (RzRw) term. However, due to the Goldstone 
s-exchange, at large quirkonium mass M the decay rate vanishes. 
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r( 1 s*+ -> vk+z) 
r( 3 s*+ -> w+z) 



Na w ot z Pwz 1 1 c w i?z 



32M 2 



Na w otz0wz 



-Rw 1 ^w — ^z J ^ (1 — Rw — Rz) 2 RwRz 
& c wRz 



64M 2 



+2{l-R w -2c 2 w R z ) 2 ( 1 



\Rs(0)\ : 



(B16) 



{l-R w ) 2 (l-R w -Rz) 2 
1 

Rz 



+ - - 4^(1 + #w + flz)) 2 W(0)|' 

tLwtlZ ) 



r( 1 p+ -> vk+z) 



RwRz 

Na w a z Pwz 
AM^{\ - R w - Rz) 2 

+4( 1- 



(l + i?^ — -Rz) 2 _|_ (l — -Rw + -Rz)' 



„2 o2 x 2 



r( 3 p+ -> w+z) 



1 — -Rw — Rz 

3 



|i?p(0)| 2 , 



i?Z 



r /3p+ ix,+ 71 = Afawgzfe [ 32c' w R 2 z 
1 1 ; IQM^l-Rw-RzfXil-Rv^ 2 



(B17) 

(B18) 
(B19) 



+ 



1 



2i?; 



8c 2 i2 z fl- 



2(\-R w ). 



1 — J?vf — Rz 



+ 



2f , ; , . Mw+Rz 

1 — -flw 



r( 3 p+ -> vk+z) = 



i?Z 

2 

|i?p(0)| 2 

iVow Ctzfiiyz 
40M 4 (l-i? w -i? z ) 



■-5— [1 — -Rw + -Rz ^ 4:CwRzY 
-Kz 



1 — i?vv — Rz 
2 

+ I 1 + 



2Rw 



2R } 



1 — -Rw — Rz ^ ^ 1 — R\v ^ 



(B20) 



r |l6[l-4r(l-ii W +iiz)] ! 



+-^-[1 + R w -Rz- 24,(1 - i? w - R z )} 2 }\R' P (0f 

K\v J 



(B21) 



4. /„/d 

Decays into two fermions only proceed via the s-channel exchange of W + . The non-zero widths with outgoing 
fcrmion masses mi _ 2 are 



16M 2 



R 2 



r(3q+,„ r J\ Na wPud 2 - Ri - R\ - R 2 + 2R X R 2 - R\ 2 
1 ( b x -*ud)- \Rs(0)\ , 



48M 2 



(l-Rw) 2 



r( 3 p+ -> ud) = 



_ NalyPud 2 - R x — R\ - R 2 + 2R X R 2 - R\ . , 



2M 4 



(1 - Rw) 2 



R> P (0)\< 



(B22) 
(B23) 
(B24) 



As expected, the 1 S partial width is proportional to m 2 /M 2 , corresponding to a chirality flip on the outgoing fermion 
line. 
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Appendix C: Decay Rates of Neutral Quirkonia 



This section summarizes the decay rates of neutral quirkonia. The decay rates differ with [19] by just a factor of 
2/3 due to a different color group. We also attempt to rewrite the decay rates so that the origins of the terms in 
the expression are manifest. In the expressions below, a i-channel quirk exchange with outgoing particles i and j 
corresponds to the factor (1 — R4 — Rj) , with i?j = m 2 /M 2 , and M is the quirkonium mass. An s-channel diagram 
exchanging particle corresponds to (1 — R^) . 



1- // 

The decays into a fermion-antifermion pair, only the s-channel 7, Z and Higgs diagram contribute. The decay of 
the 1 Sq state requires a chirality flip on the outgoing fermion line, resulting in the dependence on the fermion mass 
squared, M?, in its decay rate - similar to pseudoscalar decay. 

2 

772 

r( 1 S -> //") = 8Na%a 2 f a 2 Q p f ^r\R s (0)\ 2 , (CI) 

m z 

re* - //» - «{d + 2 «„( eae , + + («, 

rCPi -> //) = 0, (C3) 
9Nalp 3 f m 2 f , 

32Nala 2 1 /3 f 

T( 3 Pi -> //) = M4(T3flffi °ffi + (1 + 2i2/)^]|^p(0)| a , (C5) 

r( 3 P 2 //) = 0, (ce) 

where M is the quirkonium mass, az = olem / '(cjyS 2 ^), and s^ are the cosine and sine of the Weinberg angle, 
respectively, ai = and Vi = a{ — eis^ are the axial-vector and vector couplings of the Z to fermion i, with 

i = {Qi /} f° r the quirk and the outgoing ft 
velocity between the two outgoing fermions. 



* = {Qif} f° r the quirk and the outgoing fermion, respectively, Rj = m 2 /M 2 , and /3f = — 4Rf is the relative 



2. 

Only the t-channel diagram contributes decays into Zj, 



rfSo Zi) = EM M2 Z Q Q (l - Rz)\Rs(0)\ 2 , (C7) 



r( 3 ^ Z 7 ) = EA J 2 Z Q Q (l - R 2 z )\Rs(0)\ 2 , 



im 2 z 



rf^ z 7 ) = f/5 5 Q (i - i?|)|i?p(o)| 2 , 



M 2 m z 



(C8) 
(C9) 



r( 3 p h. z 7 ) = m^(i — Rz) (3 " |Ep(0)l ' (C10) 



„ &ANa.EMOtzen v n n 



(Gil) 

r( 3 p 2 -> z 7 ) = ;:;;iT/rT ( f (#1 + 3^ + e^wi 2 , (C12) 



64:Na E MO'ze'Q v Q 
5M 2 m 2 z (l - R z ) 



where the definitions of various quantities can be found in the paragraph below Eq. ( C6 ) 
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3. W+W~ 



IX 1 So 



r( 3 Si 



r( 3 p 



r( 3 p 2 



w+w-) 

W + W~) 



r( 1 Pi -> w + vk-) 
r( 3 p vf+vk") 



w+w~ 



Na 2 w P w 



8M 2 (1-2R W ) 2 

' *wPw (Rw(2- Rw) 4R W (5 + 6R W ) 



(l-2i?w) 2 
+4(1 + 20R W + 12R W ) (e Q s 2 w + - 
3Na 2 w p w 



V Q 



2Rw 
2 



V Q 



8M 2 m 2 w (l - 2R W ) 2 



Na 2 w $w 
4m^ 



1 + p% + 2R W ( 1 + 



l^5(0)| 2 , 



B 2 N 2 
Pw 



1-2R 



w 



1 - Ro 



\R' P (0)f 



1 — 2R\y 



1 — iRw + 



+2R 



w 



1 



1 — 2Ryy 



1 - 



Na 



32m 2 w 



+AR W 



J [32i& 



3-4R* 



1 — 2Rw 

(3-/3 2 w )] 2 
4 



1 — 2Rw I 1 — i?f/ 2 

i-^ Wwf, 

1 



(1 - 2i?^) 2 1 - Rz 



1 — 2R\y 1 — i?2 

2 o4 



(1 - 2R W ) 4 



\R' P (0)f 



Na 2 w l3 w 
40m^(l - 2i? w ) 2 



1 



+6P W 



2i?vi//3 4 



+8P 2 



(l-2i?w) 2 

/ B 2 , 
6+1- 



2i?vv 

1 — 2Rw J 

+ (l 



Pi 



2R\y ) 



w 

1 — 2R\y 
\R' P (0)\ 2 , 



(C13) 

(C14) 
(C15) 



(C16) 



(C17) 



(C18) 



where aw = ol em / 's 2 ^, W = VI — ^Rw is the relative velocity of the two Ws, cq,vq, and ag are the electric charge, 
vector and axial- vector couplings to the Z of the quirk, respectively, and Rw — m w /M 2 . 

4. ZZ 



rCs zz) 
r( 3 s , i -> zz) 
r( 1 Pi -> zz) 
r( 3 p -> zz) 



r( 3 Pi 



r( 3 p 2 



ZZ) 
ZZ) 



4iV( QQ + ^)a|/3| 
M 2 (1-2P Z ) 2 



|i?s(0)| ; 



3m|(l - 2i? z ) 2 

M 2 m|(l-2P Z ) 2 
2 fl _ r / 3 



l^5(0)| 2 , 
|i?p(0)| 2 , 



32m| 



32a; 



6P ? 



Q 



+8P| 



3 

1-P// 
16iVa|/3| 



32i? z v Q 
(1 - 2i?z) 2 



d-2R z y 

8(3-4R z )(a 2 Q -v 2 Q y 
(l-2Rz) 2 



)\R'f 



(0)| 2 , 



2R z v\ 

M 2 m 2 z {\-2R Z ) 2 \\-2R Z 
16Na 2 z f3 z 



-a 2 Q ) 2 |P^(0)| 2 , 



5m z 



a Q 



4 



AR% 



2 \2 



+(v 2 Q + a 2 Q ) 



4P| 



(1 - 2P Z ) 



(1-2P Z ) 2 
3 



'iRz 



2R% 



(1 - 2P Z ) 



(1-2P Z ) 2 

R' P (0)\ 2 , 



a Q 



v 2 Q 



2Rz V 

i - 2P Z y 



(C19) 
(C20) 
(C21) 

(C22) 
(C23) 

(C24) 
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where @z = \/l — ^Rz is the relative velocity between the Z's. The definitions of other quantities can be found 
below Eq. dCgb . 

5. ZP 



7Va 2 a 2 Af 2 /33 

r( 1 5 ^Zi/) = 4 ^ |i?s(0)| 2 , (C25) 

(C26) 



2Nv 2 Q all3l 
M 2 m|(l - R H - R z f 

{3 



t( i p 1 = T;;°T zg D ^ i^p(o)i 2 . (C27) 



r( d P -> ZP) = 0, (C28) 
2Na 2 n a 2 7 B ZH ( „ / R~ l 

ifPx = — / (i-P g + Pz) 2 ' z 



W T~~d~ i o 1 5 - ^ 5 i» V W(0)| a , (C29) 



1 - P z l-R H -R z 4(1 - P H - P z ) : 
3Na%a 2 7 P 5 7rT 

T ^ ~* ^ 5 M^|(l-t fi ^ ' 4(0)l - (C30) 

6. 7// 

„ Ne^OiEMOtzi^ — Rh) n 
r( 3 S x -> 7 P) = ^ V — |i?s(0)| 2 , (C31) 



6m z 



TCP, - 75) = ^^^-^V pWI 2 , 

rA Tflr, 



'Z 



HH 



Na%p H ( 9R H 6 Pl 



80to|(1 



(C32) 
(C33) 



r ( 3p 2 - pp^ - Na ^« 



^l|i?'p(0)| 2 . (C35) 
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